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Abstract 

Buckling is a significant cause of structural failure. A 
theoretical investigation has been undertaken into the dynamic 
instability of complete spherical shells which are loaded 
impulsively and made from a viscoplastic material. 

The critical mode numbers obtained in this thesis are 
similar to those obtained for cylindrical shells having the same 
R/h ratios and material parameters while the threshold velocities 
were for higher viscoplastic spherical shells than for rigid- 
plastic spherical shells. 
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Introduction 



The analysis of dynamic buckling, of elastic-plastic 
and rigid plastic cylinders has been investigated by 
Abrahamson and Goodier [1] extended by Goodier [2] , Goodier 
and Florence [3] and others [4-6]. Florence [7-8] and 
Wojewodzki [9-11] extended the analysis to viscoplastic 
cylindrical shells. V7ork on spherical shells has progressed 
through elastic-plastic and rigid plastic materials [12-14] . 

The analytical procedures of the above references were derived 
by assuming that buckling stemmed from the growth of small 
imperfections in otherwise uniform initial displacements 
and velocity fields. 

The dynamic buckling of a thin-walled complete spherical 
shell which is made from a viscoplastic material is 
investigated herein. A theoretical procedure is developed 
by first considering the unperturbed motion arising from 
a spherically symmetric initial velocity field. An 
examination is then made of the influence of small perturbations 
in this otherwise uniform velocity field and of the effect 
of small deivations from the initial sphericity of the shell. 

It transpires that certain harmonics grow very quickly and 
cause a spherical shell to exhibit a characteristic wrinkled 
shape. It is assumed that the strains which are associated 
with the perturbed state are much smaller than those 
associated with the predominant motion. Thus, no unloading 
of the shell material occurs before the time when all the 
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initial kinetic energy is absorbed by the spherical shell in 
the predominant mode. Moreover, it is further assumed that 
the wrinkled pattern in the shells is established prior to 
this time and is not modified by subsequent unloading. 
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Chapter 2 
Basic Equations 



The equilibrium equations, compatibility relations and 
strain and curvature relations have been presented in refs. [13] 
and [14] for a spherical shell. These results are repeated in 
this section for convenience. 

It may be shown for the particular case of a thin-walled 
spherical shell that the strain and curvature relations which 
were developed in reference [12] for an arbitrarily shaped 
shell, can be written 



= I " * I"'/'*' 



(la) 



eg = ^ (coseccf) V, g + cotcf) U - W + i cosec^cf) W,gW,g) (lb) 



1 

R 



e<j^g = i[cosec(|) U,g - cotcj) V + V,^ + 






(Ic) 



# (W + w,^^) 



(Id) 



• 1 • • • 

Kg = — (W + cosec^cf) W,gg + cot4> W, .) 



and K 



(j)0 



[coseccj) VJ,^Q - cot(|) coseccf) W,g] 



(j) 0 



(le) 

(l f) 



when the displacements are defined as shown in Figure 1. A 
straightforward modification of the general equations in 
ref. [12] for a shell of thickness h and density p gives the 
equilibrium equations 
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Figure 1: (a) Displacement and (b) Membrane Forces, Bending 

Moments and External Loads on a Spherical Shell. 
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- Nq) + cosec(J> N^q,0 - 

phRU + P^R =0 (2a) 

cosec4. Nq,q + + 2 cotd) N^q - 

phRV + P 2 R =0 (2b) 



and 






+ «9> + %<"'M * 

Ng cosec ^(|) W,gg + 

coseccj) N^g ■*■ ■*■ 

^^<P'(p ^'(p ^‘^'e^(pd'(p ■*■ 

^' 4 > %0'0 cosec(J> + 



COSeC^4> Ng,g 


w-e> - 


(M*. 


4><P 


2 cot 4 . 


- cotcj) 




+ 


2Mg + cosec^ 


^ ^0' 00 


+ 





2 COSeC4. M^g,^g + 

2 coseci}) coti{) M,g,g) - R^phW 

R'P3 = “ . 



+ 



(2c) 



which are consistent with equations (la) and (If) according to 
the principle of virtual work provided the influence of 
transverse shear deformation is disregarded. 

The strain and change of curvature relations (la) to (If) 
satisfy the compatibility equation 
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®e'4>4> *** ®4>'ee + ^ cot4> - 



®4>'4> + 2e^ - cosec<f> e^Q,^Q - 

• • • 

coseccf) cot(f) 6^0/0 + ■*■ ’^0^ ~ ® 



( 3 ) 



for infinitesimal displacements [12] . 

The constitutive equations for spherical shells which are 
made from rigid viscoplastic material are developed in Appendix A. 
Equations (1.8a) to (1.8f) for a spherical shell made from rigid 
viscoplastic material can be written in the form 



kh 




^4, = — (2e^ + eg) 


(4a) 


kh 




y* ^ V 


(4b) 


e .ft 

^ Y* 2 


(4c) 






kh^ 



12y’ 

kh^ 



- + Kq) 



(4d) 



and 

where 

and 



Me = 



12y’ 



kh^ 



M 






■4»e 



4)6 



12y* 2 



k = o^//3 



Y* = Y (1 - k//j7) . 



(4e) 

(4f) 

(5a) 

(5b) 



The method of solution used here is similar to that used in 
refs. [12] and [13]. It is assumed that the displacement in a 
complete spherical shell consists of dominant parts 
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(6a) 



(W, U = V = 0) and perturbation parts (w' , U' , V) so that 
W(<|>,0,T) = W(T) + W ((j>,0,T) 



and 


U(4),0,t) = U' (4), 
V(4>,0,T) = V ( 4 >, 


0,t) 

0 ,t) . 




(6b) 

(6c) 


Thus , 


the dominant strains 


and 


curvature changes are 






- — W 

® 4 ) ®0 ■ R ' 


®4>0 


= 0. 


(7a-c) 


and 


- - W 

'^ 4 ' " *^0 " ^ ' 


'^4)0 


= 0 


(7d-f) 



respectively, while the corresponding perturbation quantities 
are given by equation (la) to (If) with primes. Similarly, 



N^(<l>,0,T) = N(T) + N^((j>,0,T) (8a) 

Nq(({>,0,t) = N(T) + N^((j),0,T) (8b) 

N^0(4>,0,t) = NJq(<|),0,t) (8c) 

M^(4),0,t) = M(t) + M^((j),0,T) (8d) 

Mq(4),0,t) = M(t) + M^(4),0,t) (8e) 

and (4), 0 ,t) = ( <|) , 0 , t ) . (8f) 



If equations (6) and (8) are substituted into equation (2c) , 
then the dominant and perturbation transverse equilibrium 
equations are 

N M ^ 

2 2 — +P-- phW = 0 (9a) 

R r2 3 
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and 



R(N^ + N^) + N(w;^^ + cot<t> w;^ + cosec^ 4 ) w;^^) - 

^ ^ ^e'4. 

2 M' + cosec^<t> M',qq + 2 cosec<{) + 

2 cosec(J) cot(}) W^ 0 / 0 l ~ phR^W = 0 (9b) 

respectively, when neglecting higher order quantities. 

It may be shown that in-plane equilibrixim equations (2a) 
and ( 2 b) are satisfied identically when 



= — [F + cosec^(|) + cot(J) F,^] 






4>- 



( 10 a) 



Nfl = — IF + F, . .] 

0 j^2 <P<P 



( 10 b) 



and 



cosec4. [F,^q - cot(t> F,q] 



( 10 c) 



where F is a stress function and U=V = Pj^ = P 2 = 0. Now 
substituting the inversions of equations (4a) to (4c) and 
equations ( 12 a) to ( 12 c) and (Id) to (le) into the compatibility 
equation (3) gives 

3(V^ + 2)u,^ + [2(V^ + 1) - 1] (V^ + 2)F = 0 (11a) 



where 



u = 



W c; Fy* 

' ^ khRC ' 



C^=— T = — t 
^ p ' ^ R^ 



(llb-e) 



and 



( ) = ( ) , . . + cot(l)( ) , . + cosec^(l)( ) ,, 



(Ilf) 



(j)(j) rvy,QQ. 

The transverse equilibrixim equations (9a) and (9b) can be 
written in the form 



+ 2(1 + a^)o + P = 0 



( 12 a) 
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and 



- 3^ (V^ + 2)F - aV^u + 

a^3^[2(V2 + 1) + 1] (V^ + 2)u,^ = 0 



(12b) 



respectively, when using equations (1), (4), (6)- (8) and (10) 

. and where 



and 





kh ' 



= 



1 

12 



3 ^ 



C 

Ry* ‘ 



(12c-g) 
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Chapter 3 
Static Behavior 



Equation (12a) predicts 

P 

a (13) 

2 

when << 1, while equations (11a) and (12b) become 





[2(X„ - 1) + l)(An - 2,c^„ - 


and 


3^(X - 2)c +aX a + 

n mn n mn 

a^eM2(X^ - 1) -1U\ - 2)a^^,^ = 0 (14b) 


respectively, 


if it is assumed that 


and 


u = a P (cosd)) + 

on n ^ 

P™(cos(()) (a cos m0 + b sin m0) (15a) 

n ^ ' mn mn 

F = c^^P^ (cos({)) + 

p"*(cos({)) (c cos m0 + d sin m0) (15b) 

n ^ mn mn 


where 


X^ = n(n + 1) (15c) 



and and P^ are Legendre and associated Legendre polynomials 

of degree n and order m. Two additional equations are also 

obtained but these are identical to equations (14a) and (14b) 

provide b ^ and d„^ are substituted for a „ and If c„„ 

^ mn mn mn mn mn 

is eliminated from equations (14a) and (14b) , then 
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0 



(16a) 



3^0) a , + aX a 

n mn ' t n mn 



where 



0) = 
n 



2X - 1 

n 



+ (2X - 3) (X - 2) . 



n 



n 



(16b) 



Now, substituting equation (13) into (16a) and integrating 

Xr 



gives 



a = e 
mn 



Xn 1 c ] 

[2 in ^ J 



+ K 



(17) 



where K is a constant of integration v/hich is related to the 
initial imperfections in a shell. 3^ which equals depends 

on strain rate as y* depends on strain rate as shown in 
appendix B. With a strain rate near zero, as would be the 
case for static buckling 3^ goes to infinity but 1/3^ becomes 



Ry* R 
C C 



e y 
e ' 



+ Y 



(18) 



In the limit of very small, << y 



limit 

e ^0 
e 



R 



(19) 



where e is the initial strain rate, 
e 



becomes 



mn 



e tl 

2 o)n ^ 



Therefore equation (17) 



( 20 ) 



when e << y. 
e 



This equation is of the form of creep buckling. 
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Chapter 4 
Dynamic Buckling 



4.1 Dominant Equation 

From equation (12a) 

%'tt = <21) 

when << 1 and P = 0. 0 is not time invariant and the 

solution of the above equation will be done in two parts. 

The first part will assime the value of 0 to be constant with 
time and the second will give the exact solution, then the two 
solutions will be equated at the final displacements. For 
buckling the dominant stress is defined as 



N = - h0 



and equation (21) becomes 



u , = - 2/3 

o' TT 



( 22 ) 



(23) 



The solution of equation (23) is 




(24a) 



i« 

when satisfying the initial conditions W = 0 and W 
t = 0 and where 



■'f = 



2/3C 



V at 
o 



(24b) 



is "the time when W = 



0 . 
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The exact solution uses equation (4a) , (4b) , (7a) and (7b) 

to revise o to reflect the time dependence as follows 



kh 





(25) 



Solving for a and using equations (12c) and (2.9) yields 

n = -1- = 1 1 C 

kh “ 2 “ Ry ^o't 

-where ^ | u^,^. 

Substituting equation '(26), into (21) produces 

6C 



(26) 

(27) 



u , + — u f = — 3 

o' XT o'x 



(28) 



the equation for the exact solution. 
(28) is 



, _ 1 
^o 12 



Ry 


2 


2V 






1 + ° 




.c . 






V 



The solution of equation 
Ry 



1 - e 



6C 

___x 

Ry 



(29a) 



2C 



which satisfies the initial conditions W = 0 and W = V 

o 

and 



Ry 

T = — log 
^ 6C 



2V 



1 + 



Ry j 



at t = 0 



(29b) 



is the time when VI = 0. Setting Uq(x^) for the first solution 
equal to h^(x^) for the exact solution and solving for 
one obtains 



0^ 1 


r "N 

Ry 


2 


2V 


2V 




o _ 






— - log 


1 + ° 




0 /3 

e 


V 

oj 




Ry 


Ry 





( 30 ) 
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therefore, the dominant solution may be taken as equation (24a) 
with (24b) and (30) used with it. 

4 . 2 Perturbation Equation 

Using equations (15a-c) equation (12b) becomes 



a, -3^(1 -2)c +JLxa + 

mn TT n mn , n mn 

kh 

a232[2(X - 1) - 1] (X - 2)a , = 0 

n n mn x 



Solving equation (14a) for and substituting into 

equation (31) produces 



a , + 3^w a , + — X a = 0 

mn'xT n mn'x , , n mn 

kh 



where 



X = n(n + 1) 



n 



3(X„ - 2) 



(ij = 

n 



n 



2X - 1 



+ (2X - 



n 



3) (X - 2) 
n 



n 



= 



RY 



1 C 

+ — 



2u Ry 

ox ' 



(31) 



(32a) 

(32b) 

(32c) 



(32d) 



With a change of variable equation (32a) becomes 



a , pp - X _ 0) 
mn'^^ f n 



C C 

+ — 



UVqC ryJ 



^mn' C 



3 3V 
2 Ry 



X a = 0 
n mn 



(33a) 



where 



C = 1 - x/x. 



(33b) 
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N 

kh 



3 3V 
- + — 2 

2 RY 



(33c) 



V 



and 



u 



o' T 



(33d) 



Putting equation (33a) into a format similar to Florence [8] , 
one gets 

( n(2)l 

n ^ "r. 2 

^n a_ = 0 (34a) 



"mn'CC 



a , f. - R a =0 
mn C n mn 



where 



o'"> = S 

'n f Ry n 



(34b) 



and 



Q 



(2) = X 



n 



n f 



C 

f 2V^ n 



— + 
2 



3Vq 

Ry 



n 



(34c) 



(34d) 



As R^ is still a function of time (?) , in order to solve the 
n 

equation using Florence's solution, the value of ^ ?j '^as 

set as a constant v/ith respect to ?. A value of 3/2 was used as 
we are primarily interested near ? = 0 which negates the 
influence of the second term, thus 



R^ = 3/2 
n ' f n 



(34e) 



As given in Florence [8] the solution to equation (34a) is 



a„^(?) = A^^(?) + B^^(?)B„ 

mn mn n mn n 



( 35 ) 
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where “ 



r(2 + h^) 



exp{- (A^ + Q^) + 



1+b 



(Q K) 

n 



n 



- {(A^ + 1 + b^) + 



Q U' (f2 ) }M(fi a + { (A + 1 + b ) M(fi ) + 
n n n n n n 



M' (f2 ) } U(f2 ?) 
n n n 



(36a) 



and 



B (C) = 
mn 



r(l + a^^ + h^) 



T(2 + b ) 

' n •- 



exp{- (A^ + + 






1+b 



(fin?) 



n 






M(fi )U(Q K) 
n n 



2t, 



(36b) 



In equations (36a) and (36b), M(f2^?) and U(fij^?) are shorthand 
notations to represent Kuinmer functions [15-17] (confuent 
hypergeometric functions) of the first and second kind, 
respectively, that is 



Similarly 



M(Qj^U = M(1 + a^^ + b^, 2 + h^, Qj^C) (37a) 

U(Qj^C) = TJ(1 + a^ + bj^, 2 + h^, . (37b) 

M(Qn) = M(1 + a^ + b^^, 2 + b^, fl^^) (37c) 

U(fi)EU(l + a + b,2+b,fi) (37d) 

n n n' n' n 
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and primes over M(f2 ) and U(fi ) denote differential coefficients 

n n 

with respect to evaluated at C == 1- The remaining symbols 



are 



Q 



(2) 



a = 
n 



'n 



O 



( 1 ) 



n 



- 1 



Q 



n 



b = o(2) 



n 



n 



(37e-f) 



^ = 
n 



( 1 ) 



n 



+ 4R" 



n 



1/2 



A 



n 



0^1^ - n 

n n 



2 



(37g-h) 
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Chapter 5 

Results and Discussion 



Unfortunately, no experimental investigations appear to 
have been published on the dynamic instability of complete 
spherical shells and until reciently no theoretical work had 
been published. The results reported herein will then be 
compared with previous experimental investigations into the 
behavior of impulsively loaded cylindrical shells and in theory 
with the publications now in press. Equation (5.2) from 
Appendix E predicts mode niimbers which are 3.173/3.624 = 0.876 

times as large as Jones and Ahn's [13] estimates for similar 

h 

aluminum spherical shells and are 3.173/74'* = 1.08 times larger 
than Florence and 'Vaughan's [4] estimates for similar aluminum 
cylindrical shells. Hov/ever, as shown in Figure 2, a plot of 
the first five sets of data from Table 1, the equation for the 
critical harmonic is 

n = 2.0 (R/h) (38) 

which compares favorably with Lyons [6] curve with a slope of 
0.625. 



R/h 


Critical Harmonic n 


10 


8 


18.25 


12 


30 


15 


50 


21 


100 


32 



Table 1; A Comparison of the Critical Mode Numbers for a Visco- 
plastic Spherical Shell with Differing R/h values, 
with R = 1.460 in., Vq = 9357.6 in/sec, Y = 6500 sec 
= 46,000 psi and p = 0.0002523 Ib-sec^/in. ** . 
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Figure 2 : Critical Mode of Buckling as a Function 

of R/h V7ith = 9357.58 in/sec and 
Y = 6500 sec 
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Figure 3 shows the growth of the amplitude function with regard 
to time for the continuous spectrum of the mode numbers n. 

Showing addition figures later 7 this growth' is faster for 
larger V^, larger y and for later stages of the process 
(small C) • The change of preferred buckling mode is shown 
by a dashed line. It is evident from Table 2a that equation (38) 
provides a reasonable estimate of the preferred mode of the 
amplification function which is obtained numerically from 
the corresponding theoretical analysis using Kummer's functions 
and numerical integration of the equations. The use of the 
Kummer's functions for solution ran into problems because of 
the size of the functions and that differences of these very 
large numbers v;ere removing the first five or six significant 
figures. Even double precision was not eliminating the 
problem, for the asymptotic expansion [16] is a truncating 
series and even with double precision can't give the necessary 
significant figures. 

The Florence method gave good results that correlated 
well, up to four significant figures, with the numerical 
method. The problem with the Florence method occurred when 
the values of inputs to the Kummer function became so large 
i.e. greater than 8 that one has to use an asymptotic 
expansion. Also the last couple of computations with values 
just less than 8 one loses some accuracy. As was observed 
by Florence [8] and Wojewodzki [11] the critical modes of 
viscoplastic cylinders were constantly lower than modes of 
rigid-plastic cylinders. Likewise as shown from Table 2a, 
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Table 2a: A Comparison of the Critical Mode Numbers for a Viscoplastic Spherical Shell with the 

Corresponding Experimental Results on Cylindrical Shells. (Continued) 
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Material 


Ref. 


(ksi) 


P 

(lb~sec^/in. ** ) 


Y-l 

sec 


Steel 1015 


[1] 


102 


.0007324 


40.4 


AL 2024 T3 


[1] 


56 


. 0002523 


, 6500 


AL 6061 T6 


[1] 


44 


.0002523 


6500 


AL 6061 T6 


[5] 


44 


.0002525 


6500 


AL 6061 T6 


[6] 


46 


.0002523 


6500 


Steel 1015 


[7] 


30 


.0007324 


40.4 


Table 2b: 


values of 0^, 


p and Y 


that Correspond to 


Materials 



in Table 2a. 
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Material 

No. 


Tf 

EQ. (24b) 


C (10 ^in/sec) 
EO. (lid) 


1 


.09 4 


16.1 


2 


.110 


16.3 


3 


.108 


16.3 


4 


.085 


16.0 


5 


.158 


15.4 


6 


.140 


15.3 


7 


.182 


15.5 


8 


.146 


13.5 


9 


.153 


13.5 


10 


.166 


13.6 


11 


.173 


13.6 


12 


.171 


13.6 


13 


. 097 


13.2 


14 


.105 


13.2 


15 


.127 


13.4 


16 


.130 


13.4 


17 


.142 


13.8 


18 


.153 


13.8 


19 


.164 


13.9 


20 


.105 


13.6 


Table 2c: Values of 


and C that Correspond 


to Materials 


. in Table 2a. 


(Continued) 
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Material 

No. 


EQ. (24b) 


C (10 ^ in/sec) 
EQ. (lid) 


21 


.120 


13.7 


22 


.128 


13.7 


23 


. 200 


14.1 


24 


.179 


14.0 


25 


.169 


13.9 


26 


.179 


14.0 


27 


.123 


- 13.7 


28 


.156 


13.9 


29 


.198 


14.1 


30 


.118 


8.5 


31 


.103 


8.4 


32 


.145 


8.6 


33 


. 082 


8.3 


34 


.071 


8.2 


35 


.157 


8.7 


36 


.119 


8.5 


37 


.115 


8.5 


38 


.158 


8.7 


39 


.103 


8.4 


40 


.087 


8.3 




.076 


8.2 



Table 2c; Values of ^ and C that Correspond to Materials 
in Table 2a. 



36 



the critical modes of spherical shells are considerably lower 

than those observed experimentally on cylindrical shells with 

the same material parameters and R/h ratios which were tested 

by the authors of references [1, 5, 6, 7], with the further 

observation that correlation seems to be closer at lower R/h 

values. Of course, n, in the present case, refers to the 

★ 

critical or most amplified harmonic in the displacement 
field of a complete spherical shell which is expressed in terms 
of Legendre and associated Legendre polynomials (equation (15a) ) , 
while in reference [4] on cylindrical shells, it is associated 
with trignometric functions. 

The concept of threshold or cirtical impulse has been 
used by numerous investigators to mark the smallest impulse 
that a structure can tolerate without excessive deformation. 

It is evident from Figure 4 which is a plot of Table 3 that 
the displacement amplification function is a highly non-linear 
function of the impulsive velocity (V^) • fhe critical velocity 
is very dependent on the value of y used as shown in Figure 4 . 

It is evident from Figure 4 that the displacement 
functions at t = experience larger amplifications than do 
the velocity functions. As stated by Jones and Ahn [13] , this 
suggests that shape imperfections might exercise a much more 
significant influence on the instability of spherical shells 
than variations in the initial velocity field. 

* 

The critical harmonics are always rounded off to the nearest 
integer. 
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-4 
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Numbers in Parentheses ( ) are the Critical Modes for 



^o 

(10 in/sec) 


Y 

sec ^ 


Critical 

Mode 

.n . 


A 

mn 


Critical 
Mode 
n 


B 

mn 


5 


40 


11 


1.21 


2 


0.12 


7.5 


40 


11 


1.47 


2 


0.18 


10 


40 


11 


1.85 


5 


0.24 


11 


40 


11 


2.05 


5 


0.26 


12 


40 


11 


2.28 


6 


0.29 


15 


40 


11 


3.17 


7 


0.40 


20 


40 


12 


5.81 


8 


0.70 


30 


40 


13 


21.6 


10 


2.38 


40 


40 


14 


88.6 


12 


9.07 


45 


40 


14 


332. 


13 


28.6 


5 


6500 


11 


1.60 


2 


0.19 


7.5 


6500 


11 


2.51 


6 


0.31 


10 


6500 


12 


4.22 


8 


0.50 


11 


6500 


12 


5.28 


8 


0.61 


12 


6500 


12 


6.64 


9 


0.75 


15 


6500 


12 


13.6 


10 


1.45 


20 


6500 


13 


48.8 


11 


4.81 


25 


6500 


14 


190. 


12 


17.4 


30 


6500 


15 


788. 


13 


68.4 


32.5 


6500 


15 


1644. 


14 


138. 



Table 3; Critical Mode Numbers and Displacement Function 



Amplitudes Versus Initial Velocities at Differing 
Values of y with R = 1.460 in., R/h = 18.25, 

Oq = 46,000 psi and p = 0.0002523 Ib-sec^/in. . (Cont.) 
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V 

o 

{10 in/sec) 


Y 

-1 

sec 


Critical 

Mode 

n 


A 

mn 


Critical 

Mode 

n 


B 

mn 


5 




1,288,000 


• 11 


1.99 


5 


0.25 


7.5 




1,288,000 


12 


3.85 


8 


0.48 


10 




1,288,000 


12 


8.26 


9 


0.96 


11 




1,288,000 


12 


11.4 


10 


1.28 


12 




1,288,000 


13 


15.9 


10 


1.76 


15 




1,288,000 


13 


45.7 


11 


4.73 


17.5 




1,288,000 


14 


116. 


12 


11.4 


20 




1,288,000 


14 


304. 


13 


28.6 


25 




1,288,000 


15 


2334. 


14 


205. 


5 




900,000,000 


11 


2.24 


6 


0.29 


7.5 




900,000,000 


12 


4.89 


8 


0.61 


10 




900,000,000 


12 


11.7 


10 


1.35 


11 




900,000,000 


13 


17.5 


10 


1. 98 


12 




900,000,000 


13 


25.7 


11 


2.81 


15 




900,000,000 


14 


89.4 


12 


9.19 


17.5 




900,000,000 


14 


274. 


13 


26.6 


20 




900,000,000 


15 


856. 


13 


80.3 


22.5 




900,000,000 


15 


2902. 


14 


252. 


Table 3: 


Critical Mode Numbers and 


Displacement Function 




Amplitudes Versus 


Initial 


Velocities 


at Differing 




Values of Y with 


R = 1.460 


in. , R/h = 


^ 18.25, 




• 


0 

o 


= 46,000 psi and = 0. 


0002523 lb- 


sec /in. 


• 
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It should be noted that the present theory is an 
approximate one since the linearized from of the constitutive 
equations for rate-sensitive material vrere used. In addition 
to the various simplifications and approximations which were 
introduced in the theoretical analysis, the influence of 
transverse shear stresses and deformations have been 
disregarded. Also, the influence of unloading has been 
neglected so that the effect of possible plastic behavior 
after x = has not been considered. Moreover, the influence 
of wave disturbance has been neglected. However, with final 
times (t^) being in the neighborhood of 10-30 ysec, which can 
be calculated by using equation (lie) and Table 2c, then one 
should not neglect wave phenomena. As an example materials 
1, 10, 20, 30, 40 have values of t^ = 8.7, 13.2, 11.3, 29.1, 
and 15.9 ysec, respectively. 
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Chapter 6 
Conclusions 



A theoretical investigation has been undertaken into the 
dynamic instability of complete spherical shells which are made 
from a viscoplastic material and loaded impulsively. 

Unfortunately, no experimental results are available to examine 
the validity, or otherwise, of the various theoretical predictions. 
However, the forms of the expressions for the critical mode 
number and threshold velocity are similar to those developed 
for cylindrical shells and spherical shells, except that the 
magnitude of the numerical coefficients are different. The 
threshold velocities for viscoplastic spherical shells came out 
higher than that for rigid plastic spherical shells with the 
same R/h ratios and' material parameters. The viscoplastic theory 
predicts critical mode numbers for spherical shells which are 
similar to those observed experimentally on cylindrical shells 
with 10 < R/h < 30 but the correlation was better at the lower 
range the critical mode numbers were only one to three less 
than the rigid-plastic prediction where as the higher range 
produced mode n\ambers of four to eight less than the rigid 
plastic prediction. 



42 



Chapter 7 
Recommendations 



To further the investigation of viscoplastic effects, the 
next logical step would be to look into the non-linear 
constitutive equations. This may result in a sufficient enough 
change so as to correct some of the problems of the threshold 
velocities . 
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Appendix A 

Constitutive Equations of Viscoplastic Spherical Shells 



The constitutive equations for cylindrical shells which 
are made from a rigid viscoplastic material are developed in 



a brief derivation of constitutive equations for viscoplastic 
spherical shells. 

Consider a special linearized case of the constitutive 
equations for a rate-sensitive plastic material as derived by 
Perzyna [18] 



where e^j and are components of the plastic strain rate 

and stress tensors, k is the yield stress in simple shear and 
Y is a physical constant of the material later termed the 
viscosity constant. y is obtained by uniaxial tests as described 
later. In the above equation the function F, a function of 
strain rate and stress, is defined by F = /j^/k - 1 where 
J_ = 1/2 S. . S. . is the second invarient of the stress deviator 
with components j . In accordance with the Levy-Mises flow 
rule for rigid perfectly plastic material, as y ^2 ~ 

where k = a^//3. 

The results of uniaxial stress experiments produced the 
Cowper and Symonds [19] power-type relation 



references [8], [10], [11], and [18]. This appendix contains 



3F 



e . . = 2ky F 
ID 



( 1 . 1 ) 






( 1 . 2 ) 



Y J 
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which give the non-linearized form of function F. Inverting 
( 1 . 2 ) we get 



= Y 



- 1 



(1.3a) 



. 3 j, 

where cr and e are the equivalent stress [— S . . S . . ] and 
e e ^ 2 ID ID 

2 * * ^ 

strain rate [— e . . e. .] respectively and a is the static 
3 ID ID o 

uniaxial yield stress of the material and y and 6 are constants 
obtained experimentally for a given material 



Oe 022 ®22 ^®12 ^ ' 



• _2 r- 2 * 

^e - 7f ^^11 ^ ^ 



• 2 • 2 ^ 
11 " Z 2 ^22 ''’^12 ^ • 



(1.3b) 

(1.3c) 



Taking the linearized form of equation (1.2) we get 



a = a 
e o 



1 + 



♦ Y 

e 

e 



(1.4) 



where y is redefine as explained in Appendix E. 

In the case of plane stress ~ *^33 

equation ( 1 . 1 ) yields the follov/ing equations. 

2o. 



= 0 ) 



• Y 

11 3 






- 1 



'll ~ ^22 



(1.5a) 



22 



T 

3 






- 1 



^*^22 *^11 



(1.5b) 



^12 ^21 = 2 y 






- 1 



a 



12 






(1.5c) 
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where 



^2 = I “e' ' 



k = a //3, 
o' 



Inverting equations (1.5a-c) produces. 



11 



Y 



(2£i^ + £ 22 ) 



22 



= — ( 2£22 + 



12 



^12 ^21 



(1.5d-e) 

(1.6a) 
(1. 6b) 
(1. 6c) 



where Y* = Y(1 ~ k//j^) = y(l “ ^o^^e^ * 

Now, assume that the strain rates can be represented 



(1.6d) 
as follows: 



^11 ®1 ^ "^l 



^22 ®2 ^ *^2 



^12 ^21 ®12 ^ ’^ 12 * 



(1.7a) 

(1.7b) 

(1.7c) 



Considering that Y* is not constant and in fact is a function of 
2 , one assumes that the variation of y* with z is second order, 

a perturbation quantity. As shown in appendix B y* is dependent 

• • — .1 . ' . 

on £^ and as shown in Figure 5, £^ = £^ + £^ where £^ is a 

perturbation quantity dependent on 2 . Therefore, y* can be 

taken as a function of £ and not a function of z. Now the 

e 

constitutive equations can be integrated out as 



h/2 



-h/2 



(1.8a) 
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(1. 8b) 



h/2 



N2 = 






-h/2 

h/2 



N 



12 



^ _ kh ®12 

CJ -1 c* z ^ 

12 Y* 2 



-h/2 

h/2 



Ml = 



kh^ 



a^zdz = X--T ■ .^o 

1 i2Y* 1 2 



(2k^ + K„) 



-h/2 



(1.8c) 



(1.8d) 



h/2 



= 



02zdz 



JsiL 

12Y 



- (2k^ + K^) 



-h/2 



(1. 8e) 



h/2 



M = 
‘12 



bi2zdz 



kh^ ^12 



12Y* 2 



-h/2 



(1.8f) 
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Figure 5 ; Variation of Strain Rate Across 
Thickness of Section. 
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Appendix B 

Variations of Y and a 
o 

Equation (1.4) from Appendix A is shown in figure 6 along 
with a straight line tangent to the Cowper and Symonds 
power type relation as suggested by Florence [7] . The equation 
Florence suggested is 

“e = “Ol + "S <2.1) 

where y is the slope of the tangent line, tangent at the 

• * 

maximum initial strain rate and is the intercept of 

this tangent line at static test strain rate. The slope of 
the tangent line is obtained by differentiations of equation (1.2) 
and evaluating the results at e* 



da 



de 



• • ^ 

e = e 
e e 



j 1/6 -. 1 - 1/6 

' e 



( 2 . 2 ) 



The value of a^^ is now obtained by equation equation (2.1) and 



(1.2) evaluated at e 



• •* 



01 o 



1 + 


e 

e 


'4 


1 ] 

1 — 




. Y , 




<sj 



(2.3) 



Putting equation (2.1) into the form of equation (1.4) 



e 01 



1 + 






(2.4) 
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where now a 



01 



and Oq^/\i compare with and y from equation (1.4), 



respectively. Letting 



6 « X 

e 

• 1 e 



Y 






1/6 



1 

1 

6 



(2.5) 



and 




( 2 . 6 ) 



A replacement of and y with and respectively in the 

equations would give the identical solution with the result 
being the use of the tangent line approximation. 

The values of the material constants are 6=5 and 

Y = 40.4 sec ^ for hot rolled mild steel [20] and 6=4 and 

Y = 6500 sec ^ for 6061 T6 aluminum [21] . For thin spherical 

shells of radius Rthe initial strain rate is computed as 

• ^ 

e = 2V /R, where V is the uniform initial velocity, 
e o o . 

From Chapter 2 equation (5b) 

Y* = Y(1 - (2.7a) 

where k = 0^//3 and /j ^ (2.7b-c) 



Solving equation (1.4) for gives 



0 

o 



0 

e 




( 2 . 8 ) 



Since 1/^* is used more predominantly, substituting (2.8) into 
(2.7a) and inverting produces 




(2.9) 
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g m 
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e 
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) 



Figure 6: Cowper and Symonds Pov/er Law and Linear 

Approximation to it [11] . 
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Appendix C 

Computer Program Using Florence Method 



The following program was used to solve for ^nd 

Bmn(C)f equations (36a) and (36b), by use of Mathematics 
Library Subroutine called HYPER which is a subroutine to solve 
Kummer's function of the first kind. Problems were encountered 
when taking the difference between functions as this was 
eliminating the first five or six significant figures. 
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YMAX =0. _ FLORC032 

FL0R0033 

CONFUTE VALUES OF AN(ETA) AND BN(ETA) FLCRC034 
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COMPUTE F(l+A-G, 2-G, Z) FL0R0068 

FL0RC069 

A = 1. ♦ AT - GT FLCR0070 

G = 2. - GT FLCR0071 

Z = ZT FL0RC072 
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13 IF ( lERR.EC.O) GO TO 14 FLCR0104 

14 Rf/p = RR«=AT/GT FL0R0105 

IF(J.NE.6) GO TC 200 FLOR0106 

AMN(M) = EXP(-(ELN+0N) )*( (ELN+1.+BN)*RM0+0N*RMP)/(1.+BN) FL0R0107 

BPMP) = EXP(-(ELN + GN))«RM0/(1.+BN)*2.«TAUF FL0R0108 
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CALL NEGAM(X,GM) FLQROIAO 
G2 = GN FL0R0141 
X = AT FL0R0142 
CALL NEGAWt X,GN) FLQR0143 
G3 = GN FLGR0144 
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Appendix D 

Computer Program Using Numerical Method 



The following program solved the basic equation (34a) 

for A _(C) and B „(C) by use of Mathematics Library Subroutine 
inn inn ^ 

called DQSYK which is a subroutine to solve up to fourth order 
simultaneous differential equations v/ith initial conditions. 
This program used equation (34d) . 
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Appendix E 
Critical Impulse 



It is straightforward to show when >> 1 that A^ can 
be eliminated from equations (34c) and (34e) and using 
equation (24b) to give 



PV, 



8 



R 



n 



h/R 



/3 (a /a ) 
o e 



24/3 

^n 



- 9 



<V“e> 



1/2 



(5.1) 



where Pq/p^ is defined in equation (30). As shown in Table 4 

the values of a /a lie within the range 0.98 < a /a <1.16 

o e ■ o e 

( 2 ) 

with a medium value a /a - 1.07. Also the values of O lie 

o e 'n 

(2 ) 

within the range 1.27 < 0^ <4.42 with a medium value of 2.84. 

It is straightforward to show when using equation (34c) with 
n >> 1, that the critical harmonic 

n = 3.173(R/h)^ (5.2) 

( 2 ) 

is associated with these values of a /a and Q . 

o e n 

The threshold of dynamic plastic instability could be said 
to occur when A^j^(O) = X, where X is an arbitrary positive 
number. If X = 100, then the initial displacement imperfections 
in the critical mode are magnified by a factor of one hundred. 
This amplification is achieved when R^ - 1.7 according to 
Florence [3] asymptotic formula for equation (36a) 
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A^n<0) = e^P 



-(A„ + n ) 

n n 



(A^+l+b^)M(fi^)+n^M' (fi^) 



1 + b 



(5.3) 



n 



where the right hand quantities are defined in equations 
(37a-h) . 



(10 ^in/sec) 


(sec ) 


a /a 
o' e 


Critical 

Harmonic 


CM 

a 


5 


40 


1.00 


11 


1.27 


30 


40 


0.98 


13 


2.23 


5 


6500 


1.06 


11 


1. 35 


30 


6500 


1.03 


15 


3.99 


5 


1,288,000 


1.11 


11 


1.41 


25 


1,288,000 


1.10 


15 


4.26 


5 


900,000,000 


1.16 


11 


1.47 


22.5 


900,000,000 


1.14 


15 


4.42 



( 2 ) 

Table 4: Variation of a /a and O to various values of 

o e n 

and Y with R = 1.460 in., R/h = 18.25 = 46,000 psi 

and p = 0.0002523 Ib-sec^/in . ** . 



Finally equation (5.1) predicts a non-dimensionalized critical 
or threshold velocity 



pv. 



= 1.11 h/R 



(5.4) 



This asymptotic formula gave agreement with numerical values 
at C = 0 up to 4 significant figures. 
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This value does not compare favorably with Jones' [13] in which 
his coefficient is 17.8, since the velocities required should 
be higher than in the in rigid-plastic case. 
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